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Abstract. We discuss the relationship between the quantum Hall conductance and a fractal 
energy band structure, Hofstadter’s butterfly, on a square lattice under a magnetic field. At first, 
we calculate the Hall conductance of Hofstadter’s butterfly on the basis of the linear responce 
theory. By classifying the bands into some groups with a help of continued fraction expansion, 
we find that the conductance at the band gaps between the groups accord with the denominators 
of fractions obtained by aborting the expansion halfway. The broadening of Landau levels is 
given as an account of this correspondance. 


1. Introduction 

Through the ages, the problem of 2D Bloch electrons under an influence of vertical uniform 
magnetic field has attracted many interests of physicists. Under two limiting conditions, i.e. only 
with either magnetic field or potential, discrete Landau levels (LLs) or continuous Bloch bands 
are given as solutions BE], both of which are fundamental knowledge in solid state physics. 
A typical problem in the intermediate, or two-combined case is as follows: a single electron 
in a very strong potential under a weak magnetic field in an isotropic square lattice, where 
the Hamiltonian reduces to a proper equation named “Harper equation” [3j. While analytic 
solution does not exist for general £>, the system shows two extremely intriguing properties: 
integer quantum Hall effect (IQHE) and fractal energy band structure. The latter combines the 
adverse nature of the LLs and the Bloch bands, which is expressed in the energy diagram known 
as “Hofstadter’s butterfly” [4j. 

The IQHE is described as a quantization of Hall conductance cr, which usually occurs when 
the Fermi energy lies between the LLs. First discovered by von Klitzing et al. (5], this 
strikingly odd phenomenon has contributed to the development of condensed matter physics 
both theoretically and experimentally. For instance, since the theoretical explanation of IQHE 
was given first by Laughlin [6] based on the gauge invariance, physicists were further encouraged 
to demonstrate the validity of the method of gauge theory. Thouless et al. (or TKNN) suggested 
that a noninteracting two-dimensional electron system always shows IQHE [7] , which was later 
understood as the consequence of topological invariance of a. Namely, the Hall conductance 
can be identified as a Berry phase in the Brillouin zone (BZ) [HI 19], which is gauge invariant, 
so that a is stable unless the band gaps open and close mm- Therefore, one may expect 
that the label of each gap is related to the Hall conductance. As for the square lattice, a very 
simple yet powerful diophantic equation to associate a with the gap label was discovered by 
TKNN. However, the multivalency of the solution requires another rule to give the authentic cr, 


whose reasonable justification is still an open question. Note that, of course one may compute a 
directly from TKNN formula if needed. Other methods are Streda formula m and the adoption 
of lattice gauge theory techniques IS El. 

While the quantum Hall effect (QHE) has been a research topic, not so many works have 
focused on the fractality, an apparently interesting nature first confirmed by Hofstadter with a 
help of continued fraction expansion j3j. In this paper, we clarify the relation of the IQHE and 
the fractal structure. As a result, we find that the Hall conductance at gaps between “families” 
(the definition of a family is given in §5) can be determined without numerical calculation, whose 
uniqueness is supported by a physical background. 

The organization of this paper is as follows. §2 gives a brief derivation of Harper equation, 
which is the proper equation in our system. In §3, we estimate the conductance on the basis of 
the linear responce theory, and we describe the correspondence between the Hall conductance 
and the energy spectrum. In §4, we divide the energy spectrum into “subcells”, employing the 
notation used by Hofstadter. This enables us to state the grouping rules for families in §5. We 
discuss the relationship between two nature in §6, and the conclusion is given in §7. 

2. Harper equation 

The tight-binding Hamiltonian on an isotropic square lattice under a vertical uniform magnetic 
field is written as 


3%* = -y^Jc\cje ie ^ + (fo.c.), (1) 

{hi) 

% = ^ J* J A • dl, (2) 

where c\ (q) is a creation (annihilation) operator on the i site, and tij is a transfer integral 
between i and j sites. Peierls phase satisfies E Oij — 27t</>, with the right-hand 

plaquette 

side representing the magnetic flux per plaquette in units of magnetic flux quantum, i.e. 
B S 

cj) — — P ^/ q e uette • When we choose the Landau gauge A = (0,Hx,0), only the transfer integrals 
along the y axis acquire nonzero Peierls phase. Using the Cartesian coordinate (m, n) of site 
i (m, n=integer), the Peierls phase becomes 6ij = 2iT(f)n for the link between i = (m,n) and 
j = (m, n + 1). 

Let us assume 0 is a rational number p/q where p and q are coprime integers. In this case 
mm , the magnetic unit cell becomes q times larger in the x-direction, and correspondingly 
the BZ is q times smaller than the original one. By a straightfoward calculation, the Hamiltonian 
is transformed to 
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where each component is defined by 
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with 


c(k) = 4 (co(k),ci(k),c g _i(k)). 


(5) 
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Here, c n (k) denotes the fermion operator in the reciprocal space, 

Cm,n = C qm , +m „, n = ~^=^= ^ k ). 

V V? V L v kx , ky 

Equation |3]), which is called Harper equation, is known to have q bands when 0 = p/q. 

Figure [l] shows the energy spectrum obtained by Hofstadter [4]. It was shown m that two 



Figure 1. Energy spectrum of the tight-binding Hamiltonian of a square lattice under a 
magnetic field whose flux per plaquette is </>. The blue circle indicates a Dirac point. 


bands happen to cross linearly at zero energy, i.e. form a Dirac cone, when q is an even number. 
The blue circle shown in Fig. [l] indicates one of the Dirac cone at cf) = 1/2. Further more, the 
magnetic translation symmetry leads to a q -fold degeneracy of the Brillouin zone in the direction 
of fey, which means that each LL is q -fold degenerate Hi. In the vicinity of this Dirac point, we 
can see that the vertical uniform magnetic field measured from the Dirac point (i.e. <j) — 1/2) 
induces LLs with their energy in proportion to \Z\n\ and \/\B\ with n being the LL number. 
We refer these LLs as “the Dirac LLs”, discriminating from “the Fermi LLs” with the ordinary 
linear energy spectrum |14| . 


3. Hall Conductance 

The Hall conductance of the n-th band is computed in the linear response theory as 
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where f(E ) denotes the Fermi distribution function. Assuming zero temperature and sufficiently 
large system size, L x , L y » 1, Q is reduced to 
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where MBZ means magnetic Brillouin zone and 7 ^ is regarded as the Berry phase in the Brillouin 
zone [9]. Since the MBZ is two-torus, 7 n is integer for arbitrary n. 

As an example, we calculate the Hall conductance as a function of Fermi energy Ef at 
(f) = 9/26 as shown in Fig. [ 2 j Since there is a particle-hole symmetry, we discuss only 
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Figure 2. Hall conductance as a function of the Fermi energy for cj) — 9/26. Plateaux in family 
gaps, whose definitions are given in §5, are shown in blue lines and others in red lines. 


the region of Ef > 0 below. As Ef decreases from Ef/t — 4, the conductance changes 
as a = 0, • • • , 22,16,10,4,—2,—8,—14,—20 as shown in Fig. [ 2 | We can see that the Hall 
conductance is constant between band gaps in the energy spectrum. In the following, we clarify 
the relationship between the energy spectrum and the value of the conductance. 

4. Fractal structure of Hofstadter’s butterfly 

In order to understand the behavior of Hall conductance, we find it necessary to classify the 
magnetic Bloch bands. First, we review minimum but sufficient information of the notion 
“subcell” which was first proposed by Hofstader. We restrict our argument toO<0<l/2 since 
the identical diagram is given for the other 0 ’s. 

There are three types of subcells: L, R and C (See Fig. 3). The L and R subcells are the 
outermost bands as shown in Fig. 3(a). L n and R n (n = 0,1, 2 • • •) are in the regions of 

1 ^ 1 

-x < 0 

77 / T 3 77 / T 2 

C subcells consist of the center bands as shown in Fig. 3(a) and C n (n = 0 , 1 , 2 • • •) are in the 
regions of 

77 / , 77 / + 1 

- < <±> < - . 

2n T 1 2n T 3 

Note that the existence of gaps between the L, C and R subcells enable us to determine the 
label of each band uniquely. 
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Figure 3. (a) Butterfly divided into subcells. Red and blue lines stand for the boundary of C 
subcells and that of L, R subcells, (b) Expanded energy spectrum inside the C\ subcell. Yellow 
straight lines denote the centermost bands which have self-similar structure. 


From Fig. 3(a) we can see that each subcell, L n , i? n , C n , contains the self-similar butterfly 
pattern as in the original one for 0 < 0 < 1. Mathematically, as first claimed by Hofstadter and 
proved by MacDonald HU, an appropriate linear stretching of the local variable 0 = 1/0— [1/0] 
gives nearly the same diagram as the original butterfly pattern. Here the term ‘nearly’ reflects 
the fact that the streched diagram of L n and R n subcells differ from the original butterfly in that 
the gaps do not close at 0 = 1/2,1/4 etc, although there are no gap in the original butterfly at 
0=1/2,1/4. However, we do not need to discuss this difference further because in the following 
we restrict our argument on C subcells in which gap closes at 0 = 1/2,1/4 etc. The substructure 
inside C n can be classified again by L m , R m and C m . In the following, we denote C m in C n as 
C n C m etc; which means that u C m subcell inside C n ” Note that repetition of such operation 
determines the division of the butterfly spectrum uniquely. 

5. Grouping rule for “families” 

Next we introduce the concept of “families” in the Hofstadter’s butterfly. Figure 4 shows the 
energy spectrum obtained for 0 = 9/26. The energy levels in C\ is expanded in the r.h.s. of 
Fig. 4 and the energy bands in C\ are named in the method described in the previous section. 
We define these bands as families. 

In order to understand these families, we use T m (0) and A/(0) (m, l = integer) |4l IT : 

r ~w = 2+(m 1 +ffl - "^ <i < i2 > 

Al ^ = (7 + 2) + /T ^ 

which were introduced by Hofstadter to describe the continued fraction expansion. For example, 
0 = 9/26 can be expressed in various ways of continued fraction expansions as follows: 





















Figure 4. Energy spectrum at <j> = 
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In the general case, 0 = p/q is expressed in various ways as 
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The relation between the denominators is sin] 
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6. Relationship between conductance and energy spectrum 

We discuss the relationship between the Hall conductance and the energy spectrum on the 
basis of the ” families”. As shown in §3, the Hall conductance are, -2, -8, -14, -20 when 
0 < E f < 2. We find that the absolute values of a are equivalent to the denominators of 1/2, 
Ti(^) = 3/8, I\ 0 (^) = 5/14, ri >0>0 (i) = 7/20. In the general case, when the magnetic flux is 
4> = r mi ... mjv (/+ +1) ), we calculate r mi (|) = P 1 /Q 1 , r mim2 (|) = P 2 /Q 2 , • = 

Pn-i/Qn-i • Then the Hall conductance when Ef is located in the gap between (Z + l)th and 
Z-th family from the top is given by 


A 1 ) = \ s 9 n (fi>- -2Qi- ( 17 ) 

It is widely known that the Hall conductance is quantized in unconventional way in the gaps 

2 

of the Dirac LLs (e.g. a = t^ 4(7V + 1/2) in graphene) |19j [20]. Therefore, we can expect the 
Hall conductance to be quantized as a = 2 q(N +1/2) when infinitesimal magnetic flux A</> 
additionarily is imposed to 0 = p/q. In our case </> = 9/26, the family Ci Co Co Co corresponds to 
the n — 0 Dirac LL originated from (j) — 7/20. As antisipated from Fig. 3(b), this n — 0 Dirac 
LL is broadened as (j) decreases from 7/20, and at </> = 9/26, it becomes the family of C\ Co Co Co- 
Therefore the conductance when Ef is just above this family becomes -20 which is determined 
for (/) = 7/20. Similarly, the central three families correspond to the broadened n — 0 Dirac 
LL originated from cj) — 5/14. As a result, the conductance when Ef is just above C1C0C0L2 
becomes -14. Details would be discussed elsewhere m- 

7. Conclusion 

We have shown that the bands in 2 D electron system under a vertical uniform magnetic field 
can be classified into families, which is a notion introduced to describe the recursive structure 
of the Hofstadter’s butterfly. The Hall conductance when Ef is located in the gaps between 
families is given by the denominator of the aborted continued fraction expansion. This can be 
physically understood from the two types of LLs. 
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